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Abstract 

We perform simulations for one dimensional continuous-time random walks in two dy- 
namic random environments with fast (independent spin-flips) and slow (simple symmet- 
ric exclusion) decay of space-time correlations, respectively. We focus on the asymptotic 
speeds and the scaling limits of such random walks. We observe different behaviors de- 
pending on the dynamics of the underlying random environment and the ratio between 
the jump rate of the random walk and the one of the environment. We compare our data 
with well known results for static random environment. We observe that the non-diffusive 
regime known so far only for the static case can occur in the dynamical setup too. Such 
anomalous fluctuations give rise to a new phase diagram. Further we discuss possible 
consequences for more general static and dynamic random environments. 
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1 Introduction 



1.1 Random walk in static and dynamic random environments 

Random Walks in Random Environments (RWRE) on the integer lattice are RWs on evolv- 
ing according to random transition kernels, i.e., their transition probabilities/rates depend on 
a random field {static case) or a stochastic process (dynamic case) called Random Environ- 
ment (RE). Such models play a central role in the field of disordered systems of particles. 
The idea is to model the motion of a particle in an inhomogeneous medium. In contrast with 
standard homogeneous RW, RWRE may show several unusual phenomena as non-ballistic 
transience, non-diffusive scalings, sub-exponential decay for large deviation probabilities. All 
these features are due to impurities in the medium that produce trapping effects. Although 
they have been intensively studied by the physics and mathematics communities since the 
early 70's, except for the one-dimensional static case and few other particular situations, most 
of the results are of qualitative nature, and their behavior is far from being completely un- 
derstood. We refer the reader to \23\ [26] and [H [T3] for recent overviews of the state of the 
art in static and dynamic REs, respectively. 

In this paper we focus on two one-dimensional models in dynamic RE. In particular, the 
RW will evolve in continuous time in (two-states) REs given by two well known interacting 
particle systems: independent spin flip and simple symmetric exclusion dynamics. Several 
classical questions regarding these types of dynamical models are still open while the behavior 
of the analogous RW in a i.i.d. static case is completely understood. We perform simulations 
focusing on their long term behavior. We see how such asymptotics are influenced as a function 
of the jump rate 7 of the dynamic REs. The idea is that by tuning the speed of the REs, we 
get close to the static RE (7 close to 0) or to the averaged medium (7 approaching infinity). 
We observe different surprising phases which allow us to set some new challenging conjectures 
and open problems. 

Although the choice of the models could appear too restrictive and of limited interest, as it 
will be clear in the sequel, such particular examples present all the main features and the rich 
behavior of the general models usually considered in the RWRE literature. The conjectures 
we state can be extended to a more general setup (see Section [¥?3]l . 

The paper is organized as follows. In this section we define the models and give some 
motivation. In Section [2] we review the results known for the analogous RW in a static 
RE. Section [3] represents the main novel. We present therein the results of our simulations 
which shine a light on the behavior of the asymptotic speed (Section 13. ip and on the scaling 
limits of such processes (Sections 13.21 and 13. 3|) . When discussing each question we list several 
conjectures. In the last Section S] we present a brief description of the algorithms, we discuss 
the robustness of our numerics and possible consequences for more general RE. 

1.2 The model 

We consider a one-dimensional RW whose transition rates depend on a dynamic RE given by 
a particle system. In Section 11.2.11 we first give a rather general definition of particle systems 
and then we introduce the two explicit examples we will focus on. In Section 11.2.21 we define 
the RW in such dynamic REs. 
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1.2.1 Random Environment: particle systems 

Let 17 = {0, 1}^. Denote by Dq[0, oo) the set of paths in Q, that are right-continuous and have 
left hmits. Let {P^^r] G ^2} be a collection of probability measures on Dn[0,oo). A particle 
system 

C = (6)t>0 with 6 = {it{x) : X G Z}, (LI) 
is a Markov process on Q with law P^, when = r] gQ \s the starting configuration. Given 
a probability measure n on fi, we denote by P^{-) := /^^ P'^{-) the law of ^ when is 

drawn from /i. We say that site x is occupied by a particle (resp. vacant) at time t when 
it{x) = 1 (resp. 0). 

Informally, a particle system is a collection of particles (I's) on the integer lattice evolving 
in a Markovian way. Depending on the specific transition rates between the different config- 
urations, one obtains several types of particle systems. Each particle may interact with the 
others: the evolution of each particle is defined in terms of local transition rates that may 
depend on the state of the system in a neighborhood of the particle. For a formal construction, 
we refer the reader to Liggett [18j, Chapter L 

In the sequel we focus on two well known examples with strong and weak mixing properties, 
respectively. 

(1) Independent Spin Flip (ISF) 

Let ^ = {it)t>o be a one-dimensional independent spin- flip system, i.e., a Markov process 
on state space with generator Ljsf given by 

( W/)(??) = E [^^(^) + ^ (1 - ^(^))] [/(^") - /(^)]' ^ e (1-2) 

where A,7 > 0, / is any cylinder function on il, r/^ is the configuration obtained from r] by 
flipping the state at site x. 

In words, this process is an example of a non-interacting particle system on {0, 1}^ where 
the coordinates r]t{x) are independent two-state Markov chains, namely, at each site (indepen- 
dently with respect to the other coordinates) particles flip into holes at rate A and holes into 
particles at rate 7. This particle system has a unique ergodic measure given by the Bernoulli 
product measure with density p = "f/{"f + A) which we denote by Vp (see e.g. [IB], Chapter 
IV). 

(2) Simple Symmetric Exclusion (SSE) 

The SSE is an interacting particle system ^ in which particles perform a simple symmetric 
random walk at a certain rate 7 > with the restriction that only jumps on vacant sites are 
allowed. Formally, its generator LssE, acting on cylinder functions /, is given by 

{LssEf){v) = 7 E [/(^"'') - /(^)]' e ^' (1-3) 

where the sum runs over unordered neighboring pairs of sites in Z, and ry^'^ is the configuration 
obtained from r] by interchanging the states at sites x and y. 

It is known (see |18j , Chapter VIII) that the family of Bernoulli product measures Vp, with 
density p G (0, 1), characterize the set of equilibrium measures for this dynamics. 
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Remark 1.1. Note that the ISF and the SSE are completely different types of dynamics. 
They are both Markovian in time but while the ISF has no spatial correlations, the SSE has 
space-time correlations. The ISF has very good mixing properties due to the fact that once 7 
and X are given, no matter of the starting configuration, it will converge exponentially fast to 
the unique equilibrium given by Vp with p = 7/(7 + A). On the contrary, the SSE dynamics 
is strongly dependent on the starting configuration and therefore does not satisfy any uniform 
mixing condition. In fact, it is a conservative type of dynamics with a family of equilibria given 
by {vp : p G (0, 1)}. Because of these substantial differences, in the sequel we will informally 
say that the ISF and the SSE are examples of fast and slowly mixing dynamics, respectively. 



1.2.2 RW on particle systems 

Conditional on the particle system ^, let 

X = {Xt)t>o (1.4) 

be the continuous time random walk jumping at rate 1 with local transition probabilities 

+ 1 at rate p^t{x) + (1 - p) [1 - ^t{x)], 
x^x—1 at rate (1 - p) ^t(a;) + p [1 - 

with 

pG [1/2,1). (1.6) 

In words, the RW X jumps according to an exponential clock with rate 1, if X is on 
occupied sites (i.e. ^t{Xt) = 1), it goes to the right with probability p and to the left with 
probability 1 — p, while at vacant sites it does the opposite. 

We write Pq to denote the law of X starting from X{0) = conditional on ^, and 



Pm,o(-)=/ p^i-)p'm (1-7) 

JDn[0,oo) 

to denote the law of X averaged over 4. We refer to Pq as the quenched law and to F^^ as 
the annealed law. In what follows, when needed, we will denote by 

X{p,j,p), (1.8) 

the RW X just defined either in the ISF or in the SSE environment starting from Vp and 
jumping at rate 7. Note that in the ISF case, the parameter A is uniquely determined once 
we fix 7 and p. 

From now on we assume w.l.o.g. p G [1/2, 1). The choice of p, p G [1/2, 1) is not restrictive, 
indeed, due to symmetry, it is easy to see the following equalities in distribution 

X(p, p, 7) I X(l - p, 1 - p, 7) I -X(p, 1 - p, 7). (1.9) 



1.3 On mixing dynamics 

In our models, the REs at each site have only two possible states (0 or 1), in most of the 
literature on RWRE, the models are defined in a more general framework where infinitely 
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many states are allowed. The first paper dealing with RW in dynamic RE goes back to [9]. 
Since then, there has been intensive activity and several advances have recently been made 
showing mostly LLN, invariance principles and LDP under different assumptions on the REs 
or on the transition probabilities of the walker. See for example O El Hi El [TUl UM, [TO] (most 
of these references are in a discrete-time setting). For an extensive list of reference we refer 
the reader to [U [H] . 

One of the main difficulties in the analysis of random media arises when space-time cor- 
relations in the RE are allowed. Both models presented in Section [1.2.11 fit in this class but, 
as mentioned in Remark 11.11 their mixing properties are substantially different. The ISF 
dynamics belongs to the class of fast mixing environments which is known to be qualitatively 
similar to a homogeneous environment. In fact, a RW on this type of RE exhibits always 
diffusive scaling. 

The SSE is an example of what we called slowly mixing dynamics. For a RW driven by 
these latter types of dynamics, we are not aware of any results other than [21 [H [12]. One of 
the main result of our simulations is that the RW in (jl.4p on the SSE, similarly to the RW 
in a static RE (see Section 12. 2p , may exhibit non-diffusive behavior (see Section 13. 2p . This 
latter result is related to trapping phenomena (see Sections 12.41 and 12. 6p and suggests that, 
when considering non-uniform slowly mixing environments, the medium looks substantially 
different than a homogeneous one. This is confirmed by the rigorous annealed large deviation 
results in [2] (see the last paragraph in Section 12. 3p . Note that these results are due to the 
correlation structure of the RE but also depend on the following essential ingredients which 
produce some strong trapping effect: the one dimensional setting, the RW and the SSE being 
both nearest-neighbor, and the presence of local drifts to the right and to the left for the RW. 

1.4 Particle systems as random environments 

The reader may wonder why we consider random environments given by particle systems. 
Particle systems represent a natural physical example of a two-state dynamical RE. Particle 
system theory has been intensively developed in the last thirty years and results from this 
theory can be used in the present context (see e.g. [H [21 [Sj HI [l2l [131 [20]). Several results 
proven for such particular models can be extended to more general settings. Finally, as shown 
in this paper, these dynamics are not too complex from an algorithmic point of view allowing 
to obtain good approximations for the asymptotics. 

2 Static case and trapping phenomena: review 

We present in this section well known results for the analogous model in an i.i.d. static medium. 
Consider a static random environment rj G {0, 1}^ with law fp, the Bernoulli product measure 
with density p G [1/2, 1). Given a realization of rj, let X = {Xt)t>o be the random walk with 
transition rates (compare with (jl.Sp ) 



X — )• X + 1 at rate p7]{x) + {1 — p)[l 
X — >• X — 1 at rate {1 — p)ri{x) + p[l 



r]{x)] = c^(r?) 
rj{x)] = c~{r]) 



(2.1) 



with 



[1/2,1). 



(2.2) 
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2.1 Recurrence and LLN 



In [22] it is shown that X is recurrent when p = and transient to the right when p > ^. 
In the transient case both balhstic and non-ballistic behavior occur (see Figure [1]), namely, 
limt^oo Xt/t = ^static exists for P^, -a.e. r], and 



''^static 



= ifpG[ip], 
> if p G (p, 1]. 

In particular, for p £ (p, 1], 

I'static = Vsts,tic{p,P) = (2p - 1) 



p-p 



p{l-p) +p{l - p)' 



(2.3) 



(2.4) 




Figure 1: The sides of this square represent degenerate cases. In particular when p = 1/2 
or p = 1, the RW X does not feel anymore the environment behaving as a Simple Symmetric 
Random Walk (SSRW in the picture) or as a homogeneous RW with drift 2p — 1 (RW(p) in 
the picture), respectively. When p = 1 we are in a trivial degenerate case. When p = 1/2 we 
are in the recurrent case. Inside the square, by ()2.3p . above the diagonal we have transience 
with positive speed, while at and below the diagonal a non-ballistic transient regime holds, 
i.e. transience at zero-speed. 
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2.2 Scaling limits 



The scaling limits of one-dimensional RWRE have been derived in quite a general framework 
in a series of papers (see [T71 [2T] and ES] for a review of those results and references). 
It turns out that diffusive, super-diffusive or sub-diffusive regimes can occur. By diffusive 
regime, we mean that Xt — vt divided by y/t converges in distribution to a non-degenerate 
Gaussian distribution, while we refer to super- or sub-diffusive regime when Xt — vt has to 
be rescaled by some factor of order t", with a bigger or smaller than 1/2, respectively, to 
converge weakly to some non-degenerate distribution. We now review the different scalings in 
the i.i.d. static RE. In what follows we focus only on the annealed law. 

When X is recurrent, [21] showed that X is extremely sub-diffusive and it converges weakly 
to a non-degenerate random variable Z, namely. 



t— >oo 



(2.5) 



where cr^ is a positive constant and Z is a random variable with a non-trivial law that was 
later identified by Kesten [16j. In this case, X is called Sinai's random walk. 

When X is right-transient, |17j proved that the key quantity to determine the right scaling 
is the root s of the equation 

'c-(7?)~ 



C+(7?) 



1, 



(2.6) 



where c~{r]) and C^iri) represent the rates to jump to the left and to the right, respectively 
(see ()2.ip ). In particular, they proved that when s > 2, X is diffusive with Gaussian limiting 
distribution, while for s £ (0, 2] super- or sub-diffusivity occur with some non-trivial stable 
law of parameters {s, b) as limiting distribution (b is a constant, see also Theorem 2.3 in [26] 
for more details and references). The proof is based on the analysis of hitting times and makes 

use of the extra assumption that log 



-iv) 



has a non-arithmetic distribution. This latter is 



a delicate technical assumption (see also Remark 3 in ^]) not satisfied in our model since 



log 



■{ri) 



■iv) 



[2r/(0) - 1] log 



1 — p 
P 



(2.7) 



does have an arithmetic distribution. At the present state of the art, the role of this assump- 
tion is not entirely clear. There are examples in the literature in which by dropping it, the 
convergence in distribution does not hold (see e.g. [6] Section 8 and references therein). For 
our model ()2.ip . we performed simulations (see Figures [3] and H]) which clearly suggest that 
the arithmetic law of ()2.7p does not play any role, namely, the scaling behavior of X is like in 
the general case under the assumption of a non-arithmetic law. The following list summarizes 
the different scalings of X. 

• Diffusive: s > 2, scaling order ^/t. 

• Super-diffusive: s£ (0.5,2). 

— s G (1,2): scaling order t^^'^. 

— s = 1 : scaling order t/logt. 

— s G (0, 1): scaling order f^. 
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• Diffusive: s = 0.5, scaling order ^/i. 

• Sub-diffusive: s G (0,0.5), scaling order t^. 

In our model, the explicit solution of (|2.6p is given by 

s = s{p,p) = >0, forp,/5>l/2. (2.8) 

(¥) 

Figure [2] shows the phase diagram in {p,p) of the regimes just described. 




Figure 2: The sides of the square represent degenerate cases (see Figured]). When p = 1/2 
we are in the so called Sinai case in which X is recurrent and extremely sub-diffusive. We 
call "leaf" the region around the diagonal delimited by the curves fi{p) = (1 —pY/p^ and 
/2(p) = {p — \/p(l — p))/{2p — 1) corresponding to s = 2 and s = 0.5, respectively. The area 
above the leaf corresponds to the diffusive case with s > 2, inside the leaf we have the super- 
diffusive regime while in the lower remaining region for s G (0, 0.5) we have the sub-diffusive 
regime. 
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As we mentioned, we tested numerically the results presented so far (see Section 14.11 for 
a description of the algorithms we implemented) . Figures E] and H] show that our numerics 
match the theoretical picture just described. 




Figure 3: For each of the points marked with a black square, a dot or a cross, we performed 
numerical experiments (similar to those explained in Section \'6.'6\] to determine the scaling 
exponents of X = X{p, p). The symbols square, dot or cross mean that for the corresponding 
points, our numerical estimates gave a sub-, a super- or a diffusive scaling exponent, respec- 
tively. Note that except for a small region in between super-diffusive and diffusive regimes 
(it is reasonable to have numerical fluctuations close to a phase transition), the experiments 
confirm the theoretical picture. Figure S] provides a few explicit examples of our numerics in 
the different regimes. 
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( 1 ) Densities at s = 17.32 



( 2 ) Densities at s = 1 .64 






Figure 4: A few explicit examples of the densities obtained via our numerics in each of 
the different scaling regimes. In the right bottom picture the {p, /9)-points associated to each 
labelled plot are specified. In each plot we overlapped the densities obtained with independent 
experiments at the different times n = 2^, over samples of size M, properly rescaled. In 
particular, (1),(2),(3),(4) and (5) correspond to ballistic diffusive, ballistic super-diffusive, 
transient zero-speed super-diffusive, transient zero-speed sub-diffusive and recurrent Sinai 
case, respectively. 
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2.3 Large Deviations for the empirical speed 



For the empirical speed of a one-dimensional RW in a static RE, quenched and annealed 
Large Deviation Principles (LDPs) and refined large deviations estimates have been obtained 
in a series of papers (see e.g. [HI [151 [Ml [26]). We just mention that for the RW X defined 
through ()2.ip . when v > 0, the rate function associated to the LDP at rate t is zero on 
the whole interval [0,v]. Roughly speaking, this is saying that for 6 G [0, u), P{Xt/t w 6) 
decays sub-exponentially in t. We recall that for homogeneous RW such a decay is always 
exponential. This phenomenon is due to trapping effects which we briefly introduce in the 
next section. 

Large deviation estimates for the dynamical models in Section 11.21 were obtained in [2]. 
In particular, it is shown that in the ISF case the rate function has a unique zero (as for 
homogeneous RW) while in the SSE case, the rate function (at least under the annealed 
measure) presents a flat piece as we just described for the i.i.d. static case. 

2.4 Trapping phenomena 

The anomalous behaviors like the transient regime at zero-speed, the non-diffusivity, as well 
as the sub-exponential decay of the large deviations probabilities we reviewed, are due to the 
presence of traps in the medium, i.e., localized regions where the walk spends a long time 
with a high probability. To get an intuition, the next picture shows an explicit example of a 
trap. For a deeper insight of trapping phenomena we should introduce the random potential 
representation of the environment for which we refer the reader to the literature (see e.g. [21] 
and other references in [26]). 



Figure 5: An example of a trap, even though the global drift is to the right, i.e. p > 1/2, 
a long interval with vacant sites creates a region with local drift against the global one. To 
cross such a trap X needs an average number of trials that is exponential in the size of the 
interval. 

2.5 The averaged medium 

In the sequel we will refer to the RW in the averaged medium, y (averaged), for the homoge- 
neous nearest neighbor RW with transition rates 



where p,p G [1/2, 1). It easily follows that such a RW is right transient as soon as p > 1/2 
and p 1/2. Moreover, by the law of large numbers for i.i.d. sequences, we get that 




ooooo«««oo«« 



X — )• X -|- 1 at rate p p + {1 — p) {1 — p) 
X — >• X — 1 at rate {1 — p) p + p {1 — p) 



(2.9) 



lim 



yt( averaged) 



= V: 



t 



averae 



(2.10) 



(2p-l)(2p-l), a.s. 
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This RW would correspond to the walker defined in (jl.4p which observes at each lattice 
position a constant density of particles p. That is, between two jumps of the walker the 
environment is replaced by one that is an independent sample drawn from Vp. 

In the dynamical models of Section 11.2.11 both particle systems are assumed to be in 
equilibrium with density p and exhibit some decay of correlations in time. Hence, roughly 
speaking, if 7 approaches infinity, the environment becomes asymptotically independent be- 
tween two jumps of the walker. Therefore, under the annealed law, we expect for both models 
that there is some convergence to the averaged medium as 7 — )• 00. 

2.6 Towards the dynamic RE: dissolvence of traps 

In the previous sections we saw that the RW X in the static RE 77 G {0, 1}^ sampled from 
the Bernoulli product measure Vp presents "slow-down phenomena" due to the presence of 
traps. The dynamical models in Section 11.2.21 can be interpreted as the model in the static 
RE when we "switch on" some stochastic dynamics which allows particles of the RE to move 
(SSE dynamics) or to be created/annihilated (ISF dynamics). The natural question is then: 

How does the dynamics of the random environment influence 

the trapping effects present in the static case? 

In the sequel we will present the outcome of simulations for the empirical speed of X in 
the different REs of Section 11.2.11 and we will compare them with the static and the averaged 
medium case. Note that Vp is an equilibrium measure for both particle systems we use: ISF 
and SSE. 

At a heuristic level, one should expect that the evolution of particles in dynamic REs favors 
the dissolvence of traps, consequently the RW X in the dynamic RE should be "faster" than 
in the static medium. In other words, the long stretches of holes present at time zero, and 
responsible of the slow-down phenomena in the static case, get destroyed by the movement of 
particles in the dynamic case. As a counter effect the dynamics can create new traps during 
the evolution of the RW X. Nevertheless, in the static case, the traps are frozen, while in 
the dynamic case, all the traps have an a.s. finite survival time. Such intuitive arguments 
suggest that the displacement of the RW X should be bigger in the dynamic case than in the 
static one. Figure [6] illustrates this intuitive domination; it represents simulated trajectories 
of the RW in the three different random environments (static, ISF, SSE) starting from the 
same configuration sampled from Vp at a given p. 

Furthermore, depending on the specific dynamics of the underlying particle system, the 
survival time and the nature of a typical trap have to be different. In fact, if we consider a 
trap as in Figure [5] formed by an interval filled of holes. It is clear that in the ISF case, since 
particles can be created at each site at a given rate, such a trap gets easily destroyed. In the 
SSE cases, due to the conservation law, to dissolve such a trap, we have to wait for particles 
from outside the interval to invade it. 
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static SSE 




-50 50 1 00 1 50 200 -50 50 100 1 50 200 

location location 




-50 50 100 150 200 -50 50 100 150 200 

location location 



Figure 6: Three simulated trajectories of X in the static, ISF and SSE environments starting 
from the same configuration, {p, p) = (0.7,0.8), and 7 = 0.1 for ISF and SSE. These trajec- 
tories are compared in the bottom right picture showing the intuitive domination mentioned 
above. The background of the other plots displays realizations of the corresponding RE. In 
particular gray and white mean presence or absence of particles, respectively. Note that in 
the ISF case, due to the independence in space of the dynamics, the RE has been updated 
only around the RW trajectory, (see the description of the algorithm in Section HTT]) . 
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3 Results and conjectures 



In this section we finally present the outcome of the simulations for the asymptotics of X in 
the different cases. Section [3.11 concerns the asymptotic speed as a function of the parameters 
{PiP^l) while Sections 13.21 and 13.31 focus on the scaling limits. We give formal conjectures and 
discuss them based on the analysis of the data. The data will be presented in the form of 
figures. In the sequel, the jump rate 7 of the RE plays a central role. Throughout the paper 
we assumed, for simplicity, the RW jumping at rate 1, if instead we let it jump at any other 
rate A > 0, the same results would hold replacing 7 by 7/A. 

3.1 Asymptotic speed 

Denote by X„ = Xn{p,p,j) the position of the RW in (jl.4p after n exponential times of rate 
1. In this section we analyze the behavior of the asymptotic speed which we obtained by 
evaluating for large n 

1 *^ 

V„,=Vn{p,P,l):=^Y.^n\ (3.1) 

1=1 

over a sample of M independent experiments (the values of M's and n's will be specified in 
the figures.) 

Conjecture 3.1. Let 7 > 0, {p,p) G [1/2, 1) x [1/2, 1), and assume ^ is the SSE, then P^^ — ^ 
a.s. 

Xt{p,p,j) 
d lim =: v(jp,p,'^) G M. 

t— >oo t 

Coni ecture 13.11 should hold in great generality at least for translation invariant RE. At the 
present state of the art, the existence of an almost sure constant speed has been proven for 
dynamic REs with "good" mixing properties in space and time (see [3l O [131 [20] ) except in 
[3] which instead uses a strong elliptic condition. In particular, Conjecture 13.11 is a rigorous 
statement if i is the ISF (see [31120]'). 

Conjecture 3.2. Let ^ be either the SSE or the ISF. For any j > 0, p £ (1/2, 1), the function 
pt — > v{p,p,j) is continuous and non- decreasing. 

Note that the monotonicity is trivial once the existence of v{p,p,j) is given. Indeed, it 
follows by the fact that for any p < p', the RWs Xt{p,p,^) and Xt{p,p','y) can be coupled 
so that they remain ordered. Figure [7| below illustrates the mentioned monotonicity. In 
particular, it refers to the outcome of the numerics produced in the case of the SSE with 
p = 0.8. We remark that the same qualitative picture holds in the ISF case and for any other 
choice of p € (1/2, 1). 
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p 



Figure 7: The lower and the upper dashed curves correspond to the speeds m the static (|2.4p 
and the averaged medium ()2.10p cases, respectively. The different solid curves represent the 
function p i — > u(0.8,p, 7) at the different 7's specified. In particular, each curve is obtained 
after interpolating 11 points at distance 0.05, where each point has been obtained by averaging 
over a sample with at least M = O(IO^) simulations of the empirical speed Vn in ()3.ip with 
n = 2^^. This picture has been produced by using the SSE as RE. 

Conjecture 3.3. Let ^ he either the SSE or the ISF. Then the function 7 1 — > v{p,p,j) is 
continuous and non- decreasing. Moreover 

limv{p,p,'y) = fstatic(p,/o), (3.2) 

74-0 

lim v{p,p,-f) = i^averaged(P,/j)- (3.3) 

7— >-oo 

Conjecture 13.31 is supported by the results of the experiments presented in Figures [7] and 
[HI To state the next conjecture which describes the behavior of the speed as a function of p, 
we introduce some quantities defined in terms of 7 for fixed p G (1/2, 1). 

Conjecture 3.4. Let ^ be either the SSE or the ISF. Fix p G (1/2, 1). Define 

71 (p) := inf{7 > : v{p,p,j) > for all p > 1/2}, (3.4) 
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72(/o) := inf{7 > : p 



v{p,p,'y) is concave }, 



(3.5) 



73(/9) := inf{7 > : p i — t- v{p,p,'y) is non- decreasing }. (3.6) 
Then, the function p > — >v{p,p,^) is continuous. Moreover, 

< 7i(p) < 72(p) < 73(p) < oo. (3.7) 

Conjecture 13.41 states the existence of several critical 7's for which we see a different 
behavior of the speed as a function of p. Figure [8] below shows such a scenario in the SSE case 
at a given p = 0.8. Again, the same qualitative picture holds for any other p S (1/2, 1) or by 
considering the ISF. Above jsip) the function p 1 — > v{p,p,'y) is increasing, while it starts to 
become non- monotone for 7 < 73 (p). Below 72 (p) it looses the concavity and for 7 < 71 (p) it 
possibly starts to have a vanishing piece. A crucial issue is to understand if 

7i(/>) > 0. (3.8) 

The positivity of 71 (p) would imply a surprising transient regime with zero speed which so 
far has been proven only in the static case (see Section [2. 1|) . In the case of the ISF, it follows 
from [2Q] that 71 (p) = 0. It might still be that 71 (/o) > in the SSE case. Unfortunately, 
we feel not able to conjecture anything based on our numerics since on one hand in the 
corresponding region 7 <C 1, at the time scales we could achieve, the simulation may just be a 
weak perturbation of the static case, and on the other hand this is a very delicate phenomenon 
to test statistically since of course in any transient regime the expected speed at finite time 
is strictly positive. We therefore leave (j3.8p as a key open problem of this model. 

The loss of monotonicity for low 7's is related to the strength of the traps. In the static 
case, the speed p 1 — > v{p, p, 7) looks like the dashed lower curve in figure [5] at any fixed 
p S (1/2,1). For p big enough it starts to become decreasing until it vanishes. Intuitively, 
this is saying that when we increase p, no matter what the size of a typical trap is, the holes 
tend to act almost as reflecting barriers. 
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p 



Figure 8: p = 0.8. The lower and the upper dashed curves correspond to the speeds in the 
static (|2.4|) and the averaged medium (j2.1Up cases, respectively. The different solid curves 
represent the function p i — > v{p, 0.8, 7) at the different specified 7's. Each sohd curve is 
obtained after interpolating 26 points at distance 0.02, where each point has been obtained 
by averaging over the outcome of samples with at least M = O(IO^) independent simulations 
of the empirical speed Vn in p.ip with n > 2^^. This plot has been produced by using the 
SSE as a RE. 
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Figure [9] presents a quantitative version of Conjecture 13.41 in both the ISF and SSE cases. 
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Figure 9: Quantitative phase diagrams describing Conjecture 13.41 in the ISF and the SSE 
environments. To each (/c, 7) we associate the symbols m, c or +. They mean that the 
corresponding speed is a monotone, concave and non-concave function of p. 73 (p) would 
correspond to the curve in between the regions of m's and c's, while 72 (/o) would be the line 
separating the regions with c's and +'s. 
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3.2 Scaling limits and phase diagram: SSE 

We now turn to the analysis of the scaling hmit of X. The ISF case belongs to the class of 
dynamic RE with strong mixing properties and it is known that for such dynamics for any 
7 > 0, X satisfies a functional CLT (see e.g. [20]). In the sequel we therefore focus only on the 
SSE case. We just mention that we checked numerically the known diffusivity of the ISF case. 
We obtained excellent agreement in this case except for very few points close to degenerate 
cases, i.e. for 7 too small and (p, p) ~ (1, 1) (see also the discussion in Section |42]) . 

Our main conjecture is: 

Conjecture 3.5. Let ^ he the SSE. There exist two monotone functions 71,72 : [1/2, 1)^ 1 — > 

M+ with 

i) "fi equals zero only on the sets Ai,i = 1,2, defined by 

:= {{p, p) G [l/2,lf :p<h{p)}, (3.9) 
A2:={ip,p)G[l/2,lf :p<f2ip)}, (3.10) 
with flip) = and f2{p) = P-vgp) ^gee Figurem, 

a) 71 < 72 whenever they are non-zero, 

Hi) 7i is increasing in p and decreasing in p, for any {p,p) £ [1/2, 1)^ \ Ai,i = 1,2, 
and such that for any {p, p) S [1/2, 1)^, we have the following cases 

1. X{p,p,j) is sub-diffusive for any 7 < 7i(p, p) 

2. X(j),p,j) is super-diffusive for any 7i(p, p) < 7 < 72(p, p) 

3. X{p, p, 7) is diffusive for any 7 > 72(p, p) or at ^ = 71 (p, p) 

This conjecture is the most interesting novel result of our numerics. Figures [10] and [TT] 
show the qualitative scenario stated in Conjecture 13.51 Recall the "super-diffusive leaf" in 
Figure [21 note that the functions /i and /2 represent the lower and the upper boundary of 
the leaf, respectively. The fact that for the ISF (and more generally for RE with space-time 
correlation with exponential or fast polynomial decay) at any 7 > we have a diffusive scaling, 
can be rephrased by saying that for any 7 > the leaf vanishes and any point (p, p) G [1/2, 1)^ 
corresponds to diffusive regime. On the other hand, in case of the SSE, Conjecture 13.51 says 
that as soon as we switch on the SSE dynamics (i.e. for small 7 > 0) the leaf is still present 
and starts to move towards the p-axis as 7 increases until a certain critical 7 for which the leaf 
completely disappears. Note that the disappearance of the leaf for 7 big enough is consistent 
with the fact that as 7 increases we get closer and closer to the averaged medium case which 
is clearly diffusive. This observed phenomenon (although it could still be local, see Section 
14. 2p suggests that due to the slow-mixing properties of the exclusion dynamics, traps play a 
crucial role to determine the scaling limit of X. In particular, depending on the ratio of the 
jump rate of the walker and the one of the SSE, we can observe diffusivity or not. 



19 



Figure 10: A qualitative picture of the phase diagram described in Conjecture 13.51 in the degenerate 
cases p — 1/2 (recurrent case) and p = 1/2 (SSRW). The lower and upper dotted curves at p = 1/2 
represent 71 (p, 1/2) and 72 (p, 1/2), respectively. 




Figure 11: Two qualitative sections of the phase diagram described in Conjecture 13.51 for fixed p (on 
the left) and p (on the right). The areas of those sections in gray, black and striped correspond to 
sub-diffusive, super-diffusive and diffusive regimes, respectively. In particular, the lower and upper 
curves bounding the black areas represent 71 and 72, respectively. 
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Figures [T2HT3HT31 - [15Hl6] show some of the data supporting Conjecture 13.51 In Section 13.31 
we describe how we obtained these phase diagrams. 
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Figure 12: The section of the phase diagram in the recurrent case, i.e. p = 0.5. On the left the 
quahtative picture, on the right the outcome of our experiments supporting the qualitative phase 
diagram. The crosses, the dots, and the black squares mean that for the corresponding (p, 7) points 
our test gave an exponent equal (diffusive), bigger (super-diffusive), or smaller (sub-diffusive) than 1/2, 
respectively. More precisely, we assigned a cross at any estimated exponent in between [0.49,0.51]. 
Note that as in the lower part of the leaf for the static case (see Figure [2]), on the line dividing the 
sub- and super-diffusive regimes the scaling is diffusive. 
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Figure 13: The section of tlie piiase diagram at p = 0.7. On the left the qualitative picture, On the 
right the outcome of our experiments. As in the previous picture, crosses, dots, and squares correspond 
to diffusive, super- and sub-diffusive regimes, respectively. 
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Figure 14: The section of the phase diagram at p = 0.7. On the left the qualitative picture, on 
the right the outcome of our experiments. As in the previous pictures, crosses, dots, and squares 
correspond to diffusive, super- and sub-diffusive regimes, respectively. 
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Figure 15: Sections of data supporting the phase diagram in Conjecture 13.51 increasing in p. Crosses, 
dots, and squares correspond to diffusive, super- and sub-diffusive regimes, respectively. 
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Figure 16: Sections of data supporting the phase diagram in Conjecture 13.51 increasing in p. Crosses, 
dots, and squares correspond to diffusive, super- and sub-diffusive regimes, respectively. 
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3.3 Scaling exponents: estimating the variance for SSE 

This section is devoted to the description of the estimates we used to obtain the phase diagram 
described in Conjecture 13.51 the scahng exponents of X in the SSE environment. 

Figure [T71 below shows the same section as in Figure [12] where instead of marking "diffusive 
and non-diffusive points", we give the exphcit values of the scaling exponents. 
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Figure 17: A quantitative version of the phase diagram in Figure [T^ The numbers written at each 
point (p, 0.5,7) are the scaling exponents a*{p,p,^) that were estimated as described below. 

Each of these exponents is obtained by analyzing the variances of X at different times. 
More precisely, for each fixed triple {p,p,^), for fixed number of jumps n, we compute the 
sample standard deviation of Xn{p, p,^), namely. 



1 9 



M - 1 ^ 

i=l 



over a sample of M independent experiments, with M of order at least 10^ (the values of M 
are specified in the figures) and Vn as in (|3.ip . We then evaluate the function 

\og{SDn) 

a{n) := — 

log n 

on a log n scale, see Figures [T8HT9H20H2TH22H231 In particular, for each value of n (in the 
experiments n grows like 2^, with specified in the figures) we have performed independent 
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experiments. The proper scaling exponent is approximated by the value of a(n) with the 
biggest n which we denote by a* = a*{p,p,j) (those are the numbers in Figure [T7j) . Indeed, 
assume SDn = cn", for some positive c = c{p,p,"f), then a{n) = + a converges to a as 
n goes to infinity. 

The next figures show examples of these estimates in all described regimes. Therein we 
plotted and overlapped the empirical densities of X„ — (obtained with samples of size M) 
for different values of n to see that indeed they coincide when rescaled by . In the non- 
diffusive cases (i.e. a* ^ 1/2), we also add the plot of the same empirical densities rescaled by 
n^/^ to see that under diffusive scaling they do not coincide. In particular, they are ordered in 
time and the variances have the tendency of either vanishing (sub-diffusive) or concentrating 
(super-diffusive) . 




Figure 18: Recurrent, sub-diffusive, (p, p, 7) = (0.76,0.5,0.01). The left picture illustrates the be- 
havior of the logarithm of the standard deviation of Xn (the solid line). The dotted line represents 
the estimated exponent a* (0.76, 0.5, 0.01), while the dashed one represents the diffusive value 0.5. The 
middle plot shows that the densities obtained at different times n — 2^ rescaled by almost per- 
fectly match. In the right most plot, the same densities under diffusive scaling, looking carefully at 
the tails and around 0, one can see that they are ordered in time and tend to concentrate in as time 
increases. M denotes the sample size. The data at different times n were obtained with independent 
experiments. 
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Figure 19: Recurrent, super-diffusive, {p, p, 7) = (0.8, 0.5, 4). As in the previous figure, ttie left picture 
illustrates the behavior of the logarithm of the standard deviation of X„ (the solid line) . The dotted 
line represents the estimated exponent a* (0.8, 0.5, 4), while the dashed one represents the diffusive 
value 0.5. The middle plot shows that the densities obtained at different times n = 2^ rescaled by n" 
almost perfectly match. In the right most plot, the same densities under diffusive scaling, note that 
in this case they are ordered in time and tend to vanish as time increases. M denotes the sample size. 
The data at different times n were obtained with independent experiments. 



Log lines of standard deviations In sse: recurrent case 




binary log of numJumps 



Figure 20: Recurrent, diffusive. {p,p,j) = (0.56,0.5,0.1). As in the previous figures, the left picture 

illustrates the behavior of the logarithm of the standard deviation of X„ . The dotted line represents 
the estimated exponent a* (0.6, 0.5, 0.1). The right plot shows that the densities obtained at different 
times n = 2^ coincide under diffusive rescaling. 
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Log lines of standard deviations in sse: transient case 



Estimated scaiing: t " 0.43 



Diffusive scaling 




Figure 21: Transient, sub-diffusive. {p,p,j) = (0.74,0.55,0.001). Estimated reliable time: 
(see Section li?^. 



Log lines of standard deviations in sse: transient case Estimated scaling: t 0.57 Diffusive scaling 




Figure 22: Transient, super-diffusive, (p, p, 7) = (0.86,0.6,0.25). Estimated reliable time: 
effective running time n — 2^^ (see Section B?^ . 



28 



Log lines of standard deviations in sse: transient case 
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Figure 23: Transient, diffusive, (p, p, 7) — (0.56,0.9,0.001). Estimated reliable time: n — 2^°^ (see 
Section 112]). 

For the sake of completeness, we mention that we tried also some other methods to guess 
the right scaling exponent. A quite apt but not entirely stable method was to look for the 
exponent minimizing the total variation distance between the densities at different times. 
Another method was to fit an appropriate curve through the standard deviations directly. All 
of them gave the same qualitative picture we described here. 

4 Robustness of simulations and more general RWRE 
4.1 Algorithms 

To simulate the RW in the static RE, in the ISF and in the SSE, we implemented an algorithm 
which can be briefly described as follows. 

1) Take as INPUT (p, /), 7,n) G (0.5, 1) x [0.5, 1) x M+ x N, with n being the number of 
jumps of the RW X. 

2) Consider an interval / of size 3n centered at the starting position of the RW. 

3) Initialize the RE according to a Bernoulli product measure of parameter p. 

4) At exponential time of rate 1, update the RW position according to the underlying state 
of the RE. 

5) For each x in I, at exponential time of rate 7 update the RE at position x. 

6) Give as OUTPUT the position of the RW after n steps. 

Note that Step 5) depends on the RE. For the static case, 7 = implies that we never 
update the RE. In the ISF case, due to the independence in space of the dynamics, we update 
the RE only locally at the position of the RW X (see Figure [6]). The case of the SSE has 
been implemented using a version of the SSE either on the torus (i.e. with periodic boundary 
conditions) or by sampling at rate 7 the state of the RE at the boundaries of / from Vp (both 
approaches produce the same outcome). In particular, the size of / has been chosen of size 
3n. Particles coming from the boundaries will typically travel a distance of order ^/yn up to 
time n which guarantees that the dynamics of X is not significantly affected by those particles 
(note that the biggest 7 we considered in our numerics is of order 100). 
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4.2 Heuristics versus locality 



Conjecture 13.51 suggests a completely different scenario for the scaling limits depending on 
whether the RE is fast (ISF) or slowly mixing (SSE). This latter result is certainly the most 
interesting suggested by our numerics. Indeed, in all the papers dealing with RWDRE only 
diffusivity has been proven since, due to technical difficulties, most of the tools available are not 
suitable to treat RE presenting space-time correlations except under strong mixing conditions. 
The non-diffusivity in the experiments for the SSE could still be a local phenomenon. Namely, 
the time we run the algorithm could in principle be too small and what we are observing is only 
a perturbation of the static case or due to other local effects: e.g. presence of constants, too 
small 7, parameters in a neighborhood of the boundary of two phases or close to degenerate 
cases. Therefore, to be more confident that we are truly observing the right asymptotic 
scalings, and to have some stronger evidence in favor of Conjecture [331 we need to understand 
whether the running time of each experiment is large enough. 

We present here a heuristic argument to guess a reliable running time for small 7. The 
idea is that in each simulation the RW should perform a number of jumps sufficiently large to 
ensure that the time it takes the walker to cross a typical trap is comparable with the time 
in which this trap dissolves. 

Consider a trap of size L of consecutive holes (a stretch of size L in Z where the state 
of the RE is 0), assume the walker starts at time on the left-most hole. Let tl denote the 
number of jumps the walker needs before reaching the right-most hole of the trap. By using 
a standard gambler's ruin argument we see that 

On the other hand we can estimate the "mixing time" of the trap to be of order (7~^L)^. 
Indeed, a particle at the boundary of the trap (since it is performing a simple symmetric RW 
at rate 7) would need in average this amount of time to cross the trap. 

If E[ri] is much smaller than (j^^L)'^, this would mean that within the time the RW 
crosses the trap, the dissolvence effect due to the dynamics of the SSE is not big enough 
to play a substantial role. Therefore we would like to choose L big enough so that the two 
quantities are at least comparable, namely, we have to solve (when possible) the equation 

Once we compute (numerically) the solution L = L{p,j) of this equation, we have to be sure 
that the RW X travels a distance big enough so that it will meet with high probability such a 
big trap. The probability of observing such a disaster can be easily estimated by a geometric 
argument since the probability of a stretch of hole of size L is (1 — p)^. Therefore to be sure 
that the RW X(p, p, 7) will meet and cross at least one disaster with high probability we have 
to run the algorithm for n = n{p,p,^) steps, with n big enough so that Xn > L{1 — p)^ with 
high probability. 

It turns out that the order of such a n{p,p,^) varies a lot depending on (p,p,^). In 
most of the parameter space is unfortunately too big to be achieved in a reasonable amount 
of computing time. In the caption of Figures [^Tll22ll^ we wrote the explicit corresponding 
values of n. The fact that in most of the cases we could not run the algorithm up to such 
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an n, does not imply that the outcome of the associated experiment is only local and not 
asymptotically reliable, as the argument to deduce n is only a rough estimate. On the other 
hand, the fact that for several super-diffusive cases (as in the example in Figure [22]) we could 
achieve the associated n is a strong suggestion in favor of Conjecture 13.51 

4.3 Concluding remarks 

As mentioned before, at the current state of the art, most of the tools (regeneration, renor- 
malization, coupling, martingale approximation, etc.) developed to analyze RWRE (both 
static and dynamic REs) are still inappropriate to deal with space-time correlations unless 
the environment satisfies some uniform and fast enough mixing condition. This is not the case 
in the example of the SSE, that is why the LLN and the scaling limits for a RW driven by 
these types of REs are serious mathematical challenges. It is reasonable to think that if the 
observed unusual regimes in Conjecture 13.51 will be rigorously proven, a similar phenomenol- 
ogy can occur for other one dimensional "slowly mixing" RE (e.g. other conservative particle 
systems) providing that the RW is nearest-neighbor and has local drifts in both directions. 
In principle, we might expect some analogous scenario even in dimensions higher than one 
for static and dynamic REs presenting again some strong correlation structure. The latter 
observation is supported by the fact that a one-dimensional RW in a dynamic RE can be seen 
as a two dimensional directed RW in a static RE by interpreting the time as an extra spatial 
dimension (see e.g. [3]). 

We conclude with a table summarizing briefly some of the one-dimensional rigorous and 
numerical results we presented so far. Cells in gray mean that the corresponding result is 
non-rigorous. 





Static 


Dynamic 






SSE (slowly mixing) 


ISF (fast mixing) 


Recurrence 


p=l/2 


p = l/2 


p=l/2 


LLN 




v>0, 




with 


v>0 




v>0 






unclear for small 7 








7 big: diffusive. 




Scaling 


anomalous 


7 small: anomalous 


diffusive 


LDP rate t 


flat 


flat 


unique 


for u > 


piece 


piece 


zero 
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